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We derive a universal relation for the transverse part of triangle anomalies within 
a class of theories whose gravity dual is described by the Yang-Mills-Chern-Simons 
theory. This relation provides a set of sum rules involving the masses, decay constants 
and couplings between resonances, and leads to the formulas for the matrix elements 

■ of the vector and axial currents in the presence of the soft electromagnetic field. We 

■ also discuss that this relation is valid in real QCD at least approximately. This may 
be regarded as the anomaly matching for resonances as an analogue of that for the 

t^j- \ massless excitations in QCD. 

Oh| I. INTRODUCTION 

One distinctive feature of relativistic quantum field theories is the existence of anomalies 
[H-Ej], which is the violation of some symmetries of the classical action by quantum effects. 
In the case of global symmetries, when currents are coupled to external gauge fields, not 
00 ' an currents can be conserved. This fact is reflected in the longitudinal part of the triangle 



diagrams. The longitudinal part of triangle anomalies does not depend on the energy scale 
due to its topological nature: the triangle anomalies calculated in QCD at the level of 
quarks and ejluons are reproduced at the level of hadrons (the 't Hooft anomaly matching 
O \ condition) [4]; this leads to observable consequences for the low-energy physics involving 
pions in QCD. A well-known example is the ttq — > 2 / y decay. One can ask if the transverse 
part of the triangle graphs is also constrained. If such a constraint exists, it would have 
implications for the physics of hadron resonances (the p and a\ mesons, inparticular). 

Such a question was posed in Ref. [s| and further studied in Refs. p, 0|- It was found 
that the transverse part of the current-current correlator in an infinitesimally weak electro- 
magnetic field [denoted as wt(Q 2 ) and denned below] is not renormalized in perturbative 
QCD, and so the transverse part is related to the longitudinal part. However, chiral sym- 
metry breaking leads to a violation of this relationship. The nonperturbative aspects of the 
transverse part have been studied mostly at large Euclidean momentum Q 2 = —q 2 . Clearly, 
the main difficulty is that the transverse part of triangle anomalies has a dynamical nature 
rather than a topological one. 

In this paper, we study the transverse part of triangle anomalies using the technique of 



holography |8Hl0l] . We consider first a class of theories whose gravity dual is described by the 
Yang-Mills-Chern-Simons theory with chiral symmetry broken by boundary conditions in 
the infrared. This class of theories include the early "bottom-up" AdS / QCD model inspired 



by dimensional deconstruction and hidden local symmetry [U| and the "top-down" Sakai 



2 



Sugimoto model [12|. (Both models reproduce rather well various aspects of the physics of 
low- lying hadrons in QCD.) For models in this class, we derive the following relation for the 
transverse part of triangle anomalies: 

w T (Q 2 ) = ^ - ^ [n A (Q 2 ) - IMQ 2 )] , (i) 

for any Q 2 . Here wt{Q 2 ) is defined in Eq. ([3]) below, N c is the number of colors, f n is the pion 
decay constant, and Ha and Hy are the axial and vector current correlators, respectively. 
Equation ([T]) fully includes the nonperturbative correction and may be regarded as the 
"anomaly matching for resonances" as an analogue of that for the massless excitations 
in QCD. As will be shown, Eq. (pQ) provides a set of sum rules involving the resonance 
parameters, leading to the formulas for the matrix elements of the vector and axial currents 
in the presence of the soft electromagnetic field [see Eqs. (1541) and (|55|) ]. 

We also argue that Eq. ([T]) holds at least approximately in real QCD at both small and 
large Q 2 . 



II. TRIANGLE ANOMALIES 

First we review the triangle anomalies. We consider massless QCD with N c colors and 
Nf flavors. Let us define the correlation function of the vector current = q^^q and the 
axial current j bb = q^^^q in a weak electromagnetic background field = <9 M VJ, — d v V^ 

d a \j,jl)p = i J d*xe^(j°(x)£ b (0))p, (2) 



where t a (a = 1, 2, • • • , N 2 — 1) and t° = 1/ y2N~f are the U(iVy) flavor matrices normalized 
so that ti(t a t b ) = 5 ab /2. We also define d ab = (1/2) tr(Q{t°, t b }) where Q is the electric 
charge matrix. Since {j^Ji) P * s a Lorentz pseudo-tensor, the leading term in its expansion 
over the weak background field F^ u is a linear combination of three structures: F^ u , q^q° 'F au , 
and q v q a F aiL with F^ = (l/2)e^ l , a pF a ^ . Imposing vector current conservation q^d^il) p = 0, 
the number of independent structures reduces to two: the longitudinal and transverse parts 
with respect to q v . The general expression up to the leading order in F is 

(J»fu) P = -^ [wriq^i-q 2 ^ + q^F uv - q v q° F^) + w L (q 2 )q v q a F a ^ , (3) 

where we follow the notation of [5j. The longitudinal and transverse nature of the terms in 
this expression can be manifestly shown by using the transverse and longitudinal projection 
tensors, P^ 1 - — 77" — q^q 01 jq 2 and P"" = q^q a /q 2 '- 

(WDP = [P^T{q 2 ) + P P Jw L {q 2 )] F aj3 , (4) 



where Q 2 = —q 2 . 
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The result for wl is well-known [l|, Q] : 



w L {Q 2 ) = (5) 

This quantity does not receive corrections |3J. At the level of hadrons, the 1/Q 2 singularity 
in Eq. ([5]) is accounted for by the massless pion. 

On the other hand, the result for wt is known perturbatively j^], 

< cr \Q 2 ) = ^- (6) 



This quantity does not receive perturbative corrections as first shown by Vainshtein |5[ but 
it receives nonperturbative corrections 0, 0|. In the next section, we will show that the 
nonperturbative corrections are given in Eq. ([T|) for any Q 2 in the class of holographic QCD 
models mentioned above. 



III. HOLOGRAPHIC DESCRIPTION 



A. Setup 

The five-dimensional (5D) action of the holographic dual of our theory consists of a 
Yang-Mills (YM) and a Chern-Simons (CS) terms with a U(Afy) gauge group, 



S — Sym + Sqs 
Sym = — / d 5 xtr 



2g\z 



S, 



CS 



k / tr 



AT 2 



:A 3 T 
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A 5 



(7) 
(8) 

(9) 



Here and below, z is the fifth coordinate which runs from — zq to zq (zq > 0); the Greek 
indices /i, v, • • • denote the 4D boundary coordinates and the Latin indices M,N,- ■ • denote 
the bulk 5D coordinates. A(x, z) = Am^x m is the 5D U(iV ' f) gauge field and T = dA+iAAA 
is the field strength. They are decomposed as A = A a t a and J 7 = T a t a . 

The functions f(z) and g(z) with the conditions f(—z) = f(z) and g(—z) = g(z) (required 
by parity) are related to the metric of the bulk. For example, in the "cosh" model considered 



in [ll|, f(z) ~ cosh(,2) and g(z) = const with zq = 
f{z) ~ (1 + z 2 ) 1 ' 2 and g(z) ~ (1 + z 2 ) 1/6 with z 



oo, and in the Sakai-Sugimoto model 
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oo. In order to keep discussion general, 
we will leave f(z) and g(z) unspecified; our results below will be valid for any choice of f(z) 
and g(z) [provided that Jl° zo dz f~ 2 (z) is convergent, see Eq. ( 1231 ]. On the other hand, k will 
be fixed as k = iV c /(247r 2 ) to reproduce the correct anomaly in QCD [see Eq. ( I29p ]. In the 
top-down approach, the CS term with k = N c / (24tt 2 ) is obtained from the effective action 
of the probe D8-branes [121 ] . 
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As shown in Ref. [11] . this theory can be interpreted as a theory of mesons, which 
includes infinite towers of vector mesons and axial-vector mesons, and one massless pion. 
We decompose the gauge field A(x, z) into a parity-even part V(x, z) and a parity-odd part 

A(x,z), 

A(x,z) = V(x,z)+A(x,z), 

V(-z) = V(z), A(-z) = -A(z), 

which correspond to vector and axial- vector modes, respectively. Then boundary conditions 
are imposed at z = (which we call the IR brane): V'(0) = and A(0) = 0, where 
the derivative is taken with respect to z. Chiral symmetry is broken due to the different 
boundary conditions of V and A. The boundary conditions at z — ±zq (the UV branes) are 
the external gauge fields, 

A(z ) = A L = V + A, A(-z Q ) =A R = V~A. (11) 

Let us first recall the computation of two-point functions of currents in the absence of 
the external field F. For this purpose, the nonlinear CS term in the action can be dropped. 
We will work in the A z (x, z) = gauge. The field A^ satisfies a linear differential equation, 
which is easiest to solve in terms of the Fourier components A(x, z). The solution depends 
linearly on the boundary conditions, V® and A® , through the mode functions V(q,z), 
A(q, z), and il)(z), 

AZ(q, z) = V(q, z)P^VS Q {q) + A(q, z)P^A a a0 (q) + P^V^q) - ^(z)P^A a a0 (q) (12) 

(as will be seen later, the mode function for the longitudinal part of V is simply 1). The 
mode functions satisfy the boundary conditions 

V(q,±Zo) = l, A(q,z ) = -A(q,-z ) = l, ^{z Q ) = -ip{-z ) = 1. (13) 

The linearized field equations are given by 

d z [f(z)d z V(Q, z)] - -^-V(Q, z) = 0, (14) 

d z [f(z)d z A(Q, z)] - -^—A(Q, z) = 0, (15) 

g\z) 

d z [f(z)dMz)} = 0, (16) 

where Q 2 = —q 2 . We note that V and A are two linearly independent solutions to the same 
differential equation, so their Wronskian should be independent of z: 

f 2 (z)[V(Q,z)A'(Q,z)~A(Q,z)V'(Q,z)]=W(Q). (17) 

On the other hand, Eq. fflBT) can be solved as 

r dz' f z ° dz 

«' )=C -Lm- °-lm =1 - (18) 
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The longitundal vector mode function satisfies the same equation as Eq. (|T6j) . but with the 
boundary value of 1 at both ±z . This function is identically 1. 

Using the field equations, one can perform integration in the action by parts and the 
integral reduces to the boundary values at z — ±z : 



S 



YM 



W) 



J 2 (z)Al(q,z)d z A^(q,z) 



Z=+ZQ 



z=-z 



(19) 



Differentiating the action twice with respect to the boundary value V® , one finds the vector 
current correlation function, 



i f d 4 xe^(j;(x)j b u (0)) = 5 ab Q 2 P^U v (Q 2 



U V (Q 2 )= ±f(z)V(Q,z)V'(Q,z) 



z=-z 



—J 2 (z)V'(Q,z) 



(20) 
(21) 



and similarly for IF4. Especially, the pion decay constant f n can be obtained from the 
longitudinal part of the axial current correlation function, 



f 2 = f 2 (z)^zW(z) 



z=+zo 
\z=-zo 



2a 



or equivalent ly, 



4 

72 

J TT 



2() 



2() 



dz 



(22) 



(23) 



This expression is consistent with the one obtained in ll( as it should be. We assume that 
the right hand side of Eq. ( )23l) is convergent so that f n is finite. 



B. Longitudinal and transverse triangle anomalies 

We then take into account the effect of the CS term induced by the weak background 
field F^ u . We will work in the limit of weak background field F, and expand to linear order 
in F. For the computation of wl and wt, we can neglect the nonlinear terms in the YM 
action, because they do not include VVA interactions accompanied with e^ ua/3 tensor. 

First we note that we do not have to find the correction to the classical solution that comes 
from the CS action. Indeed, our old solution (to the Maxwell equation) is an extremum of 
the classical action, and hence a small change in the solution does not change the YM action 
to linear order. All we have to do is to substitute our old solution into the CS action. 

We then note that, unlike the YM action, the CS action is not gauge-invariant (up to 
boundaries). In order for A z = 0, we carry out the gauge transformation 

A M ->• A M - d M A, (24) 
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This keeps the transverse part of Am unchanged, but changes the longitudinal part as 
d^Az — > — d z A\. The contributions to wl and wt come from the first term in Eq. ([9]) after 
the gauge transformation: 

S cs D SKd ab F, u J d 5 x {d z V?Al - V«d z A b }). (26) 

Differentiating Sqs with respect to V® Q and A* , one obtains wl and wt- Remembering the 
definition (T4j), one has 1 

w L (Q 2 ) = J dz1/(z)V(0,z) = -^-, (27) 

MQ*) = -75- \ dzA(Q,z)V'(Q,z) } (28) 

where we took the on-shell amplitude for wi and used V(0,z) = 1. Matching between 
Eq. (127")) and the QCD result (jSJ) leads to the identification: 

K = 24^ (29) 

As seen from our derivation, wl is fixed by the boundary values alone reflecting its topological 
nature, whereas evaluating wt needs dynamical information encoded in the field equations. 
Performing the integral by parts and using Eq. f JTTj) . wt can be written as 

wt 



n n r zo 

w-wL^ VA '- AV,) 



x Nc r iz mi, ,30) 



Q 2 2Q 2 J- n / 2 W 

Using the pion decay constant (123]) . Eq. f )30|) reduces to 



On the other hand, from Eq. (I2TI) . one obtains 

2 

Q 

Combining Eqs (|3T|) and ( )32|) . one finally arrives at the relation 



U A -U V = — W(Q). (32) 



w T (Q 2 ) = ^ - ^ [n A (g 2 ) - rv(Q 2 )] , (33) 

for any Q 2 . It is clear from our derivation that this relation holds independently of f(z) 
and g(z) (i.e., the metric of the gravity). This relation for w T , which leads to the strong 



1 The expression for wl is similar to the one obtained in 13], but is different by the boundary value. 
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constraints between the resonance parameters, as we will show below, may be called the 
"anomaly matching for resonances" as an analogue of w^. 
Using the both relations for vol and wt, one also has 2 

Uffi* = -^K LR {Q 2 )PfP^F a ^ (34) 

J 7T 

for arbitrary Q 2 , where = qhl^qh is the left-handed current and j^ a = qR^^t a q R is the 
right-handed current. The form of this expression except the proportionality coefficient is 
fixed solely by the chiral symmetry SU (Nf)i, x SU(iV/)#; what we obtained here is the exact 
coefficient —N c Q 2 /(27i 2 f 2 ) including the Q 2 -dependence. 



C. Sum rules for resonances 



We shall consider the implications of the relation ( |33]) in terms of resonances (p meson, 
a\ meson, and so on). In the large N c limit, a tower of resonances with the decay widths 
T ~ 1/N C are well defined. We denote the i-th vector meson as Vi (i = 1, 2, • • - ) and j-th 
axial- vector meson as Aj (j = 1,2, •••). The wave functions for Vi and Aj in the fifth 
dimension by uA -(z) an d their masses my iA . can be determined by decomposing Eqs. (1T4"|) 
and (|T5|) into each mode with q 2 = my. Aj , respectively: 3 

(/%)'= -'fov,, W = -^» V (35) 

These functions are subject to the boundary conditions by^— z) = b v .(z), bjiA—z) = —b^iz), 
and by i (±z ) = &a,(±<zo) = with the normalization condition 

dz ——b n {z)b m {z) = 5 nm . (36) 
■zo 9 \z) 

The gauge fields V(Q,z) and A(Q,z) can be expanded as 

n<i,») = E^v,w, (37) 

A ^ z ) = E q^I b Aj (z)-^(z). (38) 



j 



Here gv t ,A are the vector and axial- vector meson decay constants defined by 

<0|#(0)|V?(p,e)) = g V r b e„ (39) 
(0|4 a (0)|A;(p,e)) = g Aj 6 ab e„ (40) 



2 In order to obtain Eq. from Eqs. ([271) and one has to add a local counter term proportional to 
a 2 F 

3 In our notation, g(z) is absorbed into bv it Aj (•?) compared with the one in |ll|: g(z)bv i , Aa , (z) — > bv i>Aj (z). 
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which can be found from Eq. (12 ID 

g Vi =- f 2 (z)b' Vi (z 



2/„W f.,\\+ z " 
20 



9Ai 



nz)b' A (z) -f 2 (z)b'( Z ) 



rZ() 



-Z() 



We also define the 7^7r-couplings g 7 v> and the 7V;A,-couplings g^Aj in 4D QCD: 

C yVi * = d ab e^g, Vi7r V«d u n b d a Vp, 
C lViAj = d ab e^g lViAj V?A b ju d a V . 

From Eq. f )26|) . these couplings are given by 4 

N c r Z0 



dzb Vi (z)ip'(z). 
| _, , dzb' v .(z)b Aj (z) 



N, r Z0 



(41) 
(42) 



(43) 
(44) 



(45) 
(46) 



2(1 



Now we are ready to write Wl and wt in terms of the resonance parameters. Substituting 
the mode expansions (J3"7j) and (138]) into Eqs. fl27j) and fl28l and performing the integration 
over z, one obtains 



47r V f 



Q 2 

4vr 2 ^ 



9Vj 

2 ' 
V, 

9Vi 



rn 



9a, 



y Q 2 + m\ Q 2 + m \. 

Therefore, Eq. ( 1271) implies the longitudinal sum rule: 

9^v i TT9v l _ N c 



(47) 
(48) 



E 



rn 



Vi 



(49) 



and Eq. (|33|) leads to the identity: 



3 Q 2 + m 2 Q 2 + m 2 A . 4tt 2 /< 



E 



#1 



m 2 (Q 2 + m 2 ) ml (Q 2 + m\ ) 



(50) 



4 Due to the identity V(0, z) = J2 k 7^ b v k 

(z) = 1, the on-shell photon in the three-point couplings can 
be replaced by the whole tower of vector mesons coupled to the photon as a manifestation of the vector 



meson dominance [l^ : 

9jVi 



9v k 

9VkVi-K—2~, 

mi, 

k v k 



9v k Viir 



J- 



dz b Vk (z)b Vz (z)^'(z), 



za 
zo 



\Ai = ^29v k v iAj ^j-, 9v k V i A j = Y^ [ dzb Vk (z)b' V[ (z)b Aj (z) 

i, m y t 47F J-zn 



k -"v k 

The quantities g^Vi-n and g^ViA, will be regarded as the "effective three-point couplings" in this respect. 
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7T 



cf = 



? 



A 6 V a 



( / 2 = m A, 



7 



(a) 



(b) 



► 



4 a 



7 



m 



7T 



7 



(C) 

FIG. 1: Diagrams contributing to the matrix elements: (a) (0| j^\^ b }p, (b) (0|j^|j4j)^,, and (c) 



for arbitrary Q 2 . Multiplying both hand sides of this identity by Q 2 + my. and then taking 
Q 2 — > —m v . limit, one obtains a set of transverse sum rules: 



E 



g-iViAjQAj 



m 



m 



,9v v , 



(51) 



for % = 1, 2, 



Similarly, 



E 



m A . - m v . 



J 9 An 



(52) 



for j 



1,2, 



J 7T 

These sum rules provide stringent constraints between the resonance 



parameters. 

These sum rules also fix the matrix elements of the vector and axial currents between 
the vacuum and one particle state (a pion, a vector meson, or an axial-vector meson) in the 
presence of the soft electromagnetic field depicted in Fig. [TJ Substituting Eqs. f )47|) and (jlSj) 
into the definitions of wl and Wt in Eq. 01]), decomposing them into the sum over i or j, 
and then using the sum rules, one finds 



(0|j"(0)|7r 6 (g))^ = iq 



v c d ab F, 



rlV 1 



(0\m\A b j (q,e)) ] 



(0\jl a (0)\V b (q,e)) } 



8 %<t 

-< i % - -V 

m A . 



Nr 



g Ai d ah F, 



2 f2 aA i 



47T 2 / 2 



— e 



-3 



'M ™2 / /| ^2^2 



m 



,/3 



9Vi - -Z7Tj-k9iV0 



4vr 2 / 2 ' 



m 



d ab F r 



a/3- 



(53) 
(54) 
(55) 
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While Eq. ( )53|) will be related to the well-known 7r — > 27 decay if one replaces the vector 
current by an on-shell photon, Eqs. (|54|) and (|55|) are the new formulas involving resonances. 
Remarkably, for fixed isospins a and b, the transverse parts of the matrix elements ( 1541) and 
( |55|) are respectively proportional to the decay constants gy i and gA j with the universal pro- 
portionality coefficient independent of species i and j (apart from the transverse projection). 
For example, for Nf = 2, one has 

AT 

3|/ 5a K>, (56) 
(0|/ 5a |^>, (7 = 1,2,...), (57) 
(0|/°|^), (i = l,2,...), (58) 



Wo! 


TV 




N 

J 7T 


(o|#W>j; 





where no summation is taken over a. We note here that the universality of the proportionality 
coefficient originates from the constant value —N c /f% with no Q 2 -dependence in front of the 
bracket in Eq. ( 1331) . 

The above sum rules and resultant matrix elements are generic to any theory with a 
Yang-Mills-Chern-Simons gravity dual in the larg e N c limit. As an example, we explicitly 



check the sum rules using the "cosh" model in Appendix |A] However, they will not 



be generally valid in a theory incorporating the scalar field corresponding to the chiral 
condensate 14H16I] (although we have a different type of sum rules which may be irrelevant 
to real QCD). We provide this counterexample in Appendix El In the next section, we will 
discuss that real QCD behaves similarly to the former class of theories with the universality 
rather than to the latter counterexample. 

If one assumes that sum rules (|5ip and fl52l) are saturated by the lowest resonances 
i — j — 1, one has 

9vi = 9a X i (59) 
9-yViAi = -4^("4-™y- (60) 



Equation (l59~j) is equivalent to the second Weinberg sum rule g v — g\ = [17| . whereas 
Eq. ( 16 Op is a new prediction. Taking experimental values for these parameters, we find 
g lp f x ~ —9.2 (and g lpa —8.0) for N c = 3. 5 This is not far from the value Ig-^/J = 7.6 ±1.1 
determined from the experimentally measured decay rate r exp (/i — > p° + 7) = 1.34 ± 0.38 



MeV 



19 



by using the formula [201 ]: 

r (/HA7) ^Kl!M^, (6i) 

24 m^rn ^ 
where d 30 = 1/4 for Q = diag(2/3, -1/3). 



5 A numerical evaluation of (|46[) using the specific metric of the Sakai-Sugimoto model gives g- IP f 1 — —3.8 
18| (after matching notation to ours), which is rather smaller than our prediction using the truncated 
sum rules. 



11 







FIG. 2: Diagrams contributing to (j^ffi, F 
respectively. 



5 v nonpcrt 



The solid and spiral lines denote quarks and gluons 



IV. REAL QCD 



Let us discuss whether the relation fl33|) is realized in real QCD. This is easy to check for 
Q 2 <C Aq CD where the dynamics is governed by the low-lying pions. Because pions do not 
contribute to %, the left hand side of f )33|) should vanish at small Q 2 . In the right hand 
side, pions only contribute to the axial correlator II a — fl/Q 2 ', the singularities of 1/Q 2 
cancel in total, and hence, Eq. is valid. 

In the opposite regime, Q 2 3> Aq CD , one can make use of the operator product expansion 
(OPE) analysis, which is an expansion of the correlator in terms of Aq CD /Q 2 . As usually 



adopted in the practical applications of the QCD sum rules [2l| . we shall neglect the a 



corrections and the anomalous dimensions of local composite operators in the OPE. Although 

our discussion 



22 



these simplifications (called the practical OPE) are numerically good 
below is approximate at this level. 

For convenience, look at the relation fl34l) instead of fl33l) . Because of the transformation 
properties under the SU(Nf)L x SXJ(Nf)n symmetry and the Lorentz symmetry, only the 
nonperturbative Lorentz pseudo-tensor condensates related to chiral symmetry breaking can 
appear in the OPE of p- The leading contributions shown in Fig. [2] read 0,0 



2 (jfiJv ) 



5 \ nonpcrt 
F 



21 

Q 6 



-q 2 O flu + q^Ovv - q u q a O, 



(62) 



where g is the QCD coupling constant and 0^ u = ((?7 At 75A a g)(g7 1/ A a g)) is the four-quark 
condensate with the SU(iV c ) color generators A a (a = 1,2, •■■ ,N 2 — 1). Using the Fierz 
transformation together with the factorization of the four-quark condensate (which can be 
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justified in the large N c limit), one has 



N 2 - 1 

CV = --^^ap(qq)(q* af3 q). (63) 



c 



24 



If we further use the magnetic susceptibility of the chiral condensate x defined by 

(?<V0) = x(qq)F„ v , (64) 

Equation fl62|) reduces to the simple form: 

= ^^^Xm 2 PfP^F aP . (65) 



On the other hand, the leading term in the OPE of IT^ is 
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a 2 N 2 - 1 o 2 

IWQ 2 ) = -^(fc7,A a g L )fe 7 ,A^)> = -^T^wY- (66) 

From Eqs. (I65p and ( I66p . that the relation (I34p holds in QCD at large Q 2 amounts to the 
condition for \ to take a special value: 

X = -^fr (67) 

Interestingly, this is the same value obtained in another way assuming the pion dominance 
in the OPE of when one turns on the quark masses 0] (see Appendix [C]) . These results 
suggest that the relation fl33l) is valid at least approximately in real QCD at both small and 
large Q 2 . 



V. CONCLUSIONS 

In this paper, we have shown a relation for the transverse part of triangle anomalies (the 
"anomaly matching for resonances") in holographic QCD. Our relation provides a set of sum 
rules involving the masses, decay constants and couplings between resonances, and leads to 
the formulas for the matrix elements of the vector and axial currents in the presence of the 
soft electromagnetic field. These results are generic to any theory with a Yang-Mills-Chern- 
Simons gravity dual where chiral symmetry is broken by the boundary conditions. 

In real QCD, our relation is also valid at least approximately when the magnetic suscep- 
tibility of the chiral condensate takes a special value x = ~ ^c/(47r 2 / 2 ). The uncertainty of 
our relation in real QCD should be resolved in the future. This is relevant to the theoreti- 
cal estimate of the hadronic electroweak contribution concerning 77*^ triangle diagrams to 
the muon anomalous magnetic moment, which can be experimentally determined to high 
precision 
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There are several open questions. Among others, it is desirable to understand our relation 
and resulting formulas for the matrix elements in the field theoretical point of view. One 



13 



can also consider its generalization to nonzero temperature and/or nonzero baryon chemical 
potential. In relation to heavy ion physics, this may lead to some possible effects on the 



"chiral magnetic effect" 
in noncentral collisions. 
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considered to explain the fluctuations of charge asymmetry 
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Appendix A: Summary of results for the "cosh" model 

In this appendix, we explicitly check our formulas in Sec. IIHI using the "cosh" model as 



an example 



11 



g{z) =95 = const, (Al) 
f(z) = -cosh(z), (A2) 

95 



and zq = oo. For completeness, we first review the results obtained in [111 ]. To match the 
notation, we assign the integer n to Vi and Aj with n = 2i — 1 for odd n and n = 2j for 
even n (due to the alternate states with the opposite parity). Then the results in 



11 



are 6 



7 / \ , ->\n -P^(tanhz) 2n + l 

= (-d-^^-J, c. = ^ ^ 7TTy , (A3) 

m 2 n = n(n + 1)A 2 , (A4) 
-A 2 



y/2n(n + l)(2n + l)—, (A5) 

95 

2A 2 

fl = % (A6) 

95 



9% 



2 " 24vr 2 ' 



(A7) 



where P^(z) are the associated Legendre functions. We then summarize the new results using 
the formulas in Sec. IIHI Introducing the variables y = tanhz and v satisfying v[y + 1) = 



6 Note that our boundary conditions (TTT|) are chosen so that the CS action is introduced in the same way 
as which are different from A(— zn) = Al and A(zq) = Ar in [11]. This entails the change of the 



sign of b n (z) (n: even) compared with 11] . 
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— Q 2 , the solutions to the field equations ([14"]) and ( !T5|) are 



K(Q,a) = -| sec(^) Vi - y 2 [^(y) + Pl{-y)i 



7T 



A(Q,*) = -secMv^^fPj^-P^-y)]. 



where sect = 1/ cost. The relations (l3Tj) and (132]) are 



n^(Q 2 )-n y (g 5 







sec 




127T 






-TrQ/A 


1 6tt 




I iv c 


A 2 


( 12tt 2 


g 2 ' 



4Q2 
A 2 

g 2 » a 2 , 
g 2 « a 2 , 



(A8) 
(A9) 



(A10) 



(All) 



w T (Q 2 



« Q 2 * L 4g 2 

~~ 71 A 2 S6C [ 2 V A 2 " 




-7TQ/A 



g 2 > a 2 , 
g 2 < a 2 , 



Therefore, the following relation is actually satisfied: 



w T (Q 2 ) = ^ - ^ [n A (g 2 ) - n v (g 5 



(A12) 



(A13) 



(A14) 



For the couplings g 7 v> and <7 7 ViA-> which we denote g inn and g 7 „ m with n = 2i — 1 and 
m = 2j, there are the "neighboring rules" : 



4^/3 

95 fir 



9-ynm = -6(n + 1} 



n(n + 2) / ( n -l)(n + l) 

:O n ,m-l + D«1 / 77; 7T77; 7TTCn,m+l- 



(A15) 
(A16) 



(2n + l)(2n + 3) ' y (2n - 1) (2n + 1) 

Using the above relations, one can easily check the longitudinal and transverse sum rules: 

g^gi _ N c _ 

N, 



2tt 2 / ot 



X ^ Q/yrvmQm 

^ m 2 m - m 5 
^ m 2 m - m 2 

n=m±l 



#5 



V2n(n + l)(2n + 1) 



#5 



y / 2m(m + l)(2m + 1) 



'4vr 2 / 2 ^' 



47r2/ 2 



(A17) 
(A18) 

(A19) 
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Appendix B: AdS/QCD with the chiral condensate 

One can test whether the relation fl33l) is realized in the AdS/QCD incorporating the 



chiral condensate 



14-16 



We consider the hard-wall model and follow the notations of 



14 



The metric is a slice of anti-de Sitter (AdS) space: 



ds 2 = — ( 

z A 



-dz 2 + dx^dXfj), 



< z < Zn 



(Bl) 



The IR cutoff z m is responsible for the confinement and fixes the scale of the p meson mass 
m p in this theory. When we are interested in the physics at large Q 2 below, we can limit 
ourselves to the region of AdS space close to the boundary and we can take the z m — > oo 
limit to simplify the computation. 

The action of the theory in the 5D bulk is 



S 
Sym 



Sym + Scs 

j5 



5, 



cs 



d x y/gtr 
K J \w b {A L ) -w 5 {A R )), 



DX\ 2 + 3\X\ 2 - -^(Fl + F j 



*9i 



(B2) 
(B3) 

(B4) 



where D^X = d^X - iA Lfl X + iXA Rfl , A L:R = A a L R r a , = d„A v - d v A^ - i[A^ A v ], and 
w(A) = AF 2 - ^A 3 F - jqA 5 . The coefficient k is fixed in Eq. (12511 . The expectation value 
of the scalar field X is determined by the classical solution as 



Xa(z) 



-m q z + -az . 



1 

2 



(B5) 



In the following, we consider the chiral limit m q = 0. 

We introduce the vector and axial-vector fields V = (A L + A R )/2 and A = (A L — A R )/2 
and we work in the V z — A z — gauge, letting V fl (q,z) = V(q, ^)l / M (g) with V being the 
source of the vector current (likewise for A^). The linearized equations of motion for the 
transverse parts V±(q, z) and A±(q,z) are 



z 



K 

z 



Q 2 

—v ± 

z 



9lv 2 



A 



0. 



0. 



(B6) 
(B7) 



with the boundary conditions V(Q, e) = A(Q, e) = 1 and V'(Q, z m ) = A'(Q, z m ) = 0. One 
can also write down the equation of motion for the longitudinal part A\\, but it is irrelevant 
to our discussion and is omitted here. 

Equation (1B6P can be solved analytically, 



V ± (Q,z)=Qz 



K l (Qz) + I l (Qz 



K (Qz m ) 

Io(Q z m) 



Z m —HX) 



>QzK x {Qz), 



(B8) 
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where K n and I n are the modified Bessel functions. Although Eq. (1B7j) does not allow for 
an analytical solution generally, one can solve perturbatively for large Q 2 , 



A ± = A + A x + . . . , 
with Aq(Q,z) = V±(Q,z). The first correction satisfies 



d 2 x A x - -d x At -A 1 = \x 4 A , 
x 



(B9) 



(BIO) 



where we define x = Qz and A = g 2 c 2 /Q 6 . The solution to this equation is given by using 
the Green's function, 



A x (x)= dx'G(x,x')\x' 4 A (x') : 



(Bll) 



where G(x,x') can be obtained from the solutions to the homogeneous part of Eq. (IBlOj) . 

= xK x (x), f 2 {x) = xh(x), (B12) 



as 



1 



G ^ x ') = ~ WTf TT7 A h{x)f2{x')0{x - x') + f 2 {x)h{x')6{x' - x)\, 
yy [Ji , h\ \ x ) 

with the Wronskian i^K^') = /1/2 — /1/2 = Using the integral, 

3 o 

dx'x' 5 Kl(x') = |, 



we find the small z behavior of Aj: 



(B13) 



(B14) 



(B15) 



This solution near the boundary is sufficient to evaluate the correlation functions below 
which are determined by the boundary values at z = e or by the integrals dominated by 
small z regions. 

The derivations of the correlation functions are similar to those in Sec. IIIII and we simply 
denote the resultant expressions here. The transverse parts of the vector and axial current 
correlation functions are 



iV(Q 2 ) 



1 Vi(Q,z) 



9lQ 2 



1 A' ± {Q,z) 



9iQ 2 



(Bi6) 

(B17) 



Since Ha(Q ) — > f%/Q f° r Q ~^ 0, the pion decay constant reads 



f 

J 7T 



1 A' ± (0,z) 

95 z 



(B18) 
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The expressions for wl and wt are 

w L {Q 2 ) = -^J* m dzA> ± (0,z)V ± (0,z) (B19) 
w T (Q 2 ) = J dzA ± (Q,z)Vi(Q,z), (B20) 

where we used Vj_(0, z) = 1. The result for wl is consistent with the anomaly matching 
condition ([5]). 7 

Now we are ready to check the validity of the relation (133]) in this theory. Let us first 
consider small Q 2 . Using V(0, z) = 1, one can easily check that Uy and wt vanish while 
TIa(Q 2 ) -»■ fllQ 2 \ thus the relation ([33]) is valid. 

On the other hand, for large Q 2 , we can expand V± and Aj_ near the boundary, 

V ± (Q,z) = l + ^(Qz) 2 \n(Q 2 z 2 ) + ..., (B21) 
A ± (Q, z) = 1 + \{Qz) 2 ln(QV) - \{Qz? g ^r + • • • , (B22) 



which lead to (up to contact terms): 



Kv(Q 2 ) = -^InQ 2 , (B23) 



where the integral 



W^-2g hQ+ 5^ (B24) 
) = ^ - — gr, (B25) 



oo 

dxx 2 K 1 (x) = 2, (B26) 

o 



is used for evaluating %. Matching the leading log behavior in Eq. (1B23P with the QCD 
result: 

U v (Q 2 ) = -^lnQ 2 , (B27) 



leads to the identification 



14 



9 12tt 2 

d = ^ (B28) 



Combining the above results, one arrives at 



MQ 2 ) = ^ - ^W^ 2 ) " U v(Q% (B29) 



7 If we take finite z m , however, the nonzero but small value of ip(z m ) at the IR brane slightly breaks the 



anomaly matching One may improve this point by adding a surface term at the IR brane 25 1. 
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for large Q 2 . Clearly, the nonperturbative correction is different from Eq. (|33|) and from the 
behavior in real QCD shown in Sec. IIV1 the coefficient in front of the bracket is Q 2 -dependent 
but not a constant —N c /f 2 . This difference originates from the OPE of wt in Eq. ( 1B25I) 
where the nonperturbative correction is proportional to 1/Q 8 rather than l/(f 2 Q 6 ). This 
will be due the absence of the field corresponding to the operator qa^q in this theory which 
is essential for the relation (133]) to be realized in real QCD at large Q 2 . One may improve 
this point by adding the tensor field H^ u corresponding to the operator qa^q in the theory, 
although it would still require a fine-tuning of parameters to reproduce the quantitatively 
correct OPE in QCD. 

In this case, one can still derive a set of transverse sum rules [but different type from 
Eqs. (l5Tj) and ( 152]) ] for highly excited resonances using Eq. (1B29j) . Since pions do not 
contribute to wt, we have only to consider the contributions from the vector and axial- 
vector mesons for wt- Similarly to Eq. (!50|) . one obtains a relation: 



; Q 2 + m 2 Q 2 + m\ 



9\ 9% 



m Vi (Q 2 + m Vi ) m % (Q 2 + 171 

W 2 , N t 



+ (B30) 



Q 2 4vr 2 ' 

for sufficiently large Q 2 . This provides a set of sum rules for highly excited states: 



£^» = 12^, «»!), 

j m Aj ~ m V t m V t 
9-yViAj9Vi - n 9A 



E^^ = 12 ^ (B31) 



They also lead to the relations for the transverse parts of the matrix elements: 

(V\jl a \V»)i = 12e a ^j-d ab F aP , (i » 1), (B32) 



9a 



(xl\j*\A% = \2e a ^-d ah F afi , (j » 1). (B33) 

These matrix elements are proportional not only to the decay constants gy i and qa^ but also 
to \jm^ /i and l/m\., respectively: there is no universality of the proportionality coefficients 
unlike Eqs. ( |5^j) and ( 1551) . This again comes from the 1/Q 2 behavior in front of the bracket in 
Eq. (1B29|) . and is different from real QCD where this factor should be approximately replaced 
by a constant value. Therefore, we expect that real QCD would have the properties 
and dSSJ) rather than (1B321) and flB33|) . 



Appendix C: Magnetic susceptibility of the chiral condensate 

In this appendix, we review the derivation of the magnetic susceptibility of the chiral 
condensate by Vainshtein 0. Let us consider the modifications of the longitudinal part of 
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the correlator (J4]) when we turn on the degenerate quark masses m q . For Q 2 3> Aq CD , the 
leading contribution can be found using the OPE, 



N c ppa _2m q (qa^q) 



Q 2 



47T 



pa± p0\\p 



2 /J f 



a/3 



4tt 2 q 2 ^ q 4 

2iV c levr^rn, (gg) 



Q 2 



Q 4 



+ o 



i 



(Cl) 



where we used the definition of x i n Eq. flBTj) . For Q 2 <C Aq CD , the pion propagator is 
replaced by the massive one: 



Q n ,,± 



4vr 2 



paj-pp\\ p 



2N r 



pa±p/3\\p 



Q 2 + ml 
~2N C 2N c m 2 n 



Q 2 



Q 1 



(C2) 



If one assumes the extrapolation of the l/Q 4 term in the bracket in Eq. ( 1C2I) to large Q 2 to 
be matched against that in Eq. (1C1I) . one finds 



X 



N c ml 



lir 2 m g (qq}' 

Using the Gell-Mann-Oakes-Renner relation for pions, 

fl m l = -2m,<gg), 

Equation (1C3I) reduces to 



X 



(C3) 



(C4) 



(C5) 
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